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ABSTRACT 

In the standard model of cosmology, structure emerges out of non-rotational 
flow and the angular momentum of collapsing halos is induced by tidal torques. 
The growth of halo angular momentum in the linear and quasi-linear phases is 
associated with a shear, curl-free, flow and it is well described within the linear 
framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as 
haloes approach turn around and virialization. At that stage the flow field around 
halos has non-zero vorticity. Using a cosmological simulation, we have examined 
the importance of the curl of the velocity field (vorticity) in determining halo 
spin, finding a strong alignment between the two. We have also examined the 
alignment of vorticity with the principle axes of the shear tensor, finding that it 
tends to be perpendicular to the axis along which material is collapsing fastest 
(ei). This behavior is independent of halo masses and cosmic web environment. 
Our results agree with previous findings on the tendency of halo spin to be 
perpendicular to e 1; and of the spin of (simulated) halos and (observed) galaxies 
to be aligned with the large-scale structure. Our results imply that angular 
momentum growth proceeds in two distinct phases. In the first phase angular 
momentum emerges out of a shear, curl-free, potential flow, as described by 
TTT. In the second phase, in which haloes approach virialization, the angular 
momentum emerges out of a vortical flow and halo spin becomes strongly aligned 
with the vorticity of the ambient flow field. 

Subject headings: Dark Matter: haloes 
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Introduction 



Most galaxies, especially late type spirals, are observed to be rotating. Any robust 
theory of structure and galaxy formation needs to thus account for the origin of angular 
momentum. In physics in general and hydrodynamics in p articular, rotation is often as- 
sociat ed with vorticity. Not surprisingly, the early ideas of iHeisenberg fc von Weizsacker 
(119481 ) attempted to attribute the rotation of spiral nebulae to a primeval vorticity. The 
so-called theory of cosmic turbulence suggested that density perturbations, out of which 
structure evolves, emerged out of pri mordial turbulence and th at galaxies got their rotation 
from a prim ordial vortical flow (e.g. lOzernoi fc Cherninl Il968l . cf. the excellent review of 
Joneslll976l ). These ideas were superseded by an alternative theory which has survived and 
emerged as the current standard model of cosmology: structure grows by the evolution of 
smal l primordial density perturbations via gravita tional instability in an expanding universe 



(e.g. iPress fc Schechterlll974l ; IWhite &: Reeslll978l ). By construction the primordial flow field 
is a potential flow (thus curl-free) and therefore an alternative model fo r the origin of the 
gala ctic sp i n need s to be invoked. The first su ch ideas were d eveloped by IStrombergi (119341 ) 
and iHoyld (119511 ) and then independently by iPeeblesI (Il969l ). who suggested that tidal in- 
teraction between collapsing haloes may be responsible for galaxy spins. A ccordingly, th e 
induced internal velocity field in the torquing objects, is a curl-free shear flow (lBinneylll974l ). 
It follows that the resulting angular momentum of cosmological objects is associated with a 
shear and not a vortical flow. 

The classi cal picture of the origin of angular momentum was expanded by Tidal Torque 
Theory (TTT, IWhitdll984f ). This seminal paper forms the base for the standard lore of the 
origin of angular momentum in modern cosmology. TTT provides a framewor k for calculating 



the dynamics and statist i cs of the growth of the angula r momentum (e.g. iHoffmanl 11986 
Heavens fc Peacock! Il988l ; ISteinmetz fc Bartelmannl|l995l ). In fact, it is the only framework 
that enables an analytical analysis of the growth of the angular momentum of dark matter 
(DM) halos. It does so by considering the growth of the Lagrangian mass of DM halos 
within the linear theory of gravitational instability. Yet TTT is of limited validity when 
consid ering highly non-linear virialized halos. Using numerical simulations iPorciani et al. 



( 120021 ) examined the growth of the angular momentum of the Lagrangian mass associated 
with z = halos and found that TTT provides a useful description until around ~ 60% of the 
effective turn-around time. At z = the relative error in the predictions for the magnitude 
of the spin is of the order of unity and the mean scatter in the misalignment between the 
simulated and the TTT predicted direction of the spin is ~ 50°. It is clear that non-linear 
effects weaken the applicability of TTT on the scale of virial halos. 



A close examination of TTT reveals that one of the model's most important tenets is the 
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fact that the primordial flow field is curl-free and therefore that angular momentum emerges 
out of a shear flow. This principle clearly do es not hold deep i n the non-linear regime, 
wher e the cosmic flow fi eld becomes vortical (IKnebe et al.l 120061 ; iPueblas &: Scoccimarro 
20091 ; iKitaura et al.l 120121 ) . This has prompted us to examine the possibility that the spin 
of halos is a result of the vorticity of the cosmic flow field on scales greater than the virial 
radius. This is line with recent efforts to understand the emergence of the cosmic web and 
the p roperties of DM halos in terms of the local differential properties of the cosmic flow 
field Jrloffman et al.lboij ILibeskind et al.ll2012aB ) 



Connecting galaxy or halo properties to their envi ronments is no t a new idea, and has 



its roots in the seminal density-morphology relation o f iDresslerl (I198C ). A number of recent 



papers [jAragon-Calvo et al.ll2007l ; iBrunino et al.ll2007t iHoffman et al 



2012: Libeskind et al. 



2012a| jbl. among others) have investigated the relation between halo spin and shape with the 
LSS in numerical simulations using a variety of different techniques to define the cosmic web. 



Libeskind et al.l (j2012al Jbl) used the velocity shear tensor to show that halo spin tends 



to be perpendicular to the eigenvector corresponding to the largest eigenvalue of the shear - 
the direction of greatest collapse (or, in voids, of weakest expansion). It is therefore natural 
to extended these ideas to the anti-symmetric part and examine how the direction of the 
curl of the velocity field is distributed relative to the eigenvectors of the shear tensor and 
halo spin. 



Methods 



The same DM only simulation is used here as in IHoffman et al.l ( 120121 ) . It follows 
the evolution of 1024 3 particles in a 64/i -1 Mpc periodic box employing standard WMAP5, 
ACDM cosmological parameters with a mass resolution of ~ 1.89 x lO^ ^Mm and a spatial 
resolution of lft. -1 kpc. The halo finder AHF (IKnollmann &: Knebell2009l ) is run on the z = 
particle distribution. Around ~ 10 6 haloes more massive than lO lo /i -1 M are considered 
here. 



2.1. Vorticity and shear calculation 

The velocity field is evaluated using the Clouds-In-Cell (CIC) approach calculated on a 
256 3 grid giving a spatial resolution of 250/i -1 kpc. The CIC gridding of the velocity field has 
been performed after excising the DM particles that make up virial halos. The rationale for 
that is to avoid 'double counting' effects when relating the properties of halos with the their 
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ambient velocity field. Differentiation of the velocity field is performed in Fourier space, using 
an FFT. The CIC velocity field is Gaussian smoothed with kernel lengths of r g = 250/i~ 1 kpc, 
500/i -1 kpc and lOOO/i^kpc. A given smoothing length is then matched to the appropriate 
halo mass been, such that r g roughly equals four virial radii. 

The velocity field v at each point in space may by expanded as the sum of symmetric 
and anti-symmetric components. The symmetric component is the shear tensor: 

1 (dv a dv@\ 



2H \dr/3 dr c „ 
x, y, and z and H is the Hubble constant. Since E c 



(1) 



where a, f3 — x, y, and z and H is the Hubble constant. Since T, a/3 is symmetric its eigen- 
values (Ai, A2, and A3) are real and its trace is equal to the divergence of the velocity, which 
is proportional to the overdensity in the linear regime. The number of positive (negative) 
eigenvalues corresponds directly to the number of axes along which matter is contracting 
(expanding). If zero, one, two or three axes are simulta neously collapsing, th e region may 



be classified as a void, sheet, filament or knot (e.g. see iHoffman et al.l 120121 ) . These axes 
are just the eigenvectors of the shear tensor (e 1; e 2 and e 3 ). The anti-symmetric component 
of the expansion is the vorticity oj defined as the curl of t he velocity fie l d: cu = V x v. 
Vorticity generation is a strongly non-linear effect, shown by iKitaura et al.l (120121 ) to appear 
only in third order Lagrangian perturbation theory. Initial perturbations are assumed to 
be effectively curl-free since in an expanding un iverse conservation o f angular momentum 
implies that any vortical flow will decay linearly (IBouchet et al.lll995l ). 



Haloes are divided into three mass bins: low mass (10 10 - 10 11 h" l M Q ), intermediate mass 
(10 11 — lO 12 /i~ 1 M ) and high mass (10 12 — \Q 1?, h~ l M Q ). The median halo virial radius (and 
standard deviation) of each mass bin are 54±13/i _1 kpc, 115±30/i _1 kpc, and 252±64/i _1 kpc, 
respectively. The alignment between halo spin and the vorticity is examined on scales > 4r v i r 
using the smallest (250/i _1 kpc), intermediate (500/i _1 kpc) and largest (1000/i _1 kpc) gaussian 
smoothing kernel for the smallest, intermediate and largest mass bin, respectively. 



3. Results 

3.1. Alignment of cosmic vorticity with the velocity shear. 

In this section the alignment of vorticity with the shear tensor eigenvectors is examined. 
The probability distribution of the angle formed between the vorticity vector and each of the 
three eigenvector of the shear tensor at each grid point in the simulation is shown in FigUfa). 
Since the eigenvectors of the shear denote axes (not directions), the interval is confined to 
I cos fi\ G [0, 1]. Fig[T]^a) clearly shows that the vorticity tends to strongly align perpendicular 
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to ei and in the e 2 — e3 plane, evidenced by the highly non uniform distributions. The 
apparent parallel alignment with e 2 and e 3 is driven by the perpendicularity with respect to 
ej. Note that there is a slight tendency for the vorticity to be more closely aligned with e3, 
than e 2 , indicating a non- uniform azimuthal distribution. 

The significance of the alignment is quantified by calculating the average offset between 
the measured distribution and a uniform one, in units of the Poisson error a per bin of a 
uniform distribution. If less than unity, the measured distribution is consistent with random. 
All three alignments in Fig. Ufa) show a ~ 6.5a statistical significance indicating they are 
fully inconsistent with a random distribution. The median angles between the vorticity and 
ei, e 2 and e 3 are 0.29, 0.60 and 0.54, respectively. The vorticity of 'the velocity field tends to 
be perpendicular to the main axis of the shear and lie in the plane of the intermediate and 
minor axis. 

The alignment between the vorticity and the eigenvectors of the shear tensor is also 
shown for each cosmic web environment as the thin lines in Fig. Ufa). The perpendicular 
alignment between u and ei decreases in strength from sheets, to filaments, to knots. This 
is likely due to the increased erratic behavior of these vectors deep in non-linear regimes. 
Regardless of which web element is examined, u tends to cluster more towards e 3 than e 2 
and always away from ei. 

Given the dominance of the shear tensor in shaping the local dynamics it is only natural 
to study the orientation of the vorticity with respect to its three eigenvectors. Since the 
eigenvectors are non-directional lines (axes) the coordinate system they define does not span 
the full three-dimensional volume, but just the positive octant of the cartesian grid. Each 
vorticity vector can thus be normalized to unity, and its direction projected onto a unit 
sphere defined by the eigenvectors. 

The distributions of u on (one eighth of) the unit sphere can be seen in an equal 
area projection in Fig. |2j The perpendicularity with respect to ei is clearly visible as the 
vorticity tends to inhabit regions close to the e 2 -e3 plane (the "equator"). As exemplified 
by the stronger parallel alignment between u and e 3 in Fig. [TTa), at low latitude (i.e. for u's 
that are far from ej), the vorticity is closer to e 3 than to e 2 . At higher latitudes, the locus 
of vorticity vectors is more uniform moving closer to e 2 at the highest latitudes. 

In Fig. [3] we present the vorticity and shear field around a generic overdensity in a 
10 x 10/t _1 Mpc region. In overdense regions (where the vorticity is well defined) ei tends to 
point in the same direction across many Mpc, indicating that the axis of greatest collapse is 
correlated on these scales. In the same region, the perpendicularity of the vorticity and ei 
is clearly visible. 
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3.2. Alignment of cosmic vorticity with halo spins. 



Since iLibeskind et al.l ( l2012af ) showed that halo spin axes tend to align with the prin- 
ciple axes of the shear tensor, and in the preceding section a correlation between the shear 
eigenvectors and the vorticity was shown, the alignment of vorticity and halo spin is now 
examined. Such an alignment would be indicative of a rotational flow, rather than the clas- 
sical picture that haloes form out of a potential, irrotational flow. In Fig. H^b) we show the 
alignment of halo spin with the vorticity for three halo mass bins. Note that in this case the 
angle is between two vectors and as such allowed to range from [-1,1]. 

Fig. 1(b) shows a clear correlation of halo spin with vorticity for all mass bins. The 
strength of this alignment is mass dependent, with larger halo spin being more tightly 
aligned with the vorticity: the smallest, intermediate and highest mass bins have medi- 
ans of cos /j, =0.11, 0.22, and 0.37. For the smallest mass bin, the haloes have the weakest 
alignment, yet display the strongest signal (~ 5.3a) statistically speaking, since this mass bin 
contains the largest number of haloes. The highest mass bin is inconsistent with a random 
distribution at the ~98% confidence level and shows the strongest alignment. Halo spin is 
strongly correlated with the large scale vorticity of the velocity field. 



4. Conclusions and discussion 

Using a high resolution cosmological simulation the relationship between halo spin and 
the relevant directions that characterize the large scale structure (LSS) is studied. The 
velocity field has been evaluated on a smoothed, clouds-in-cells grid. The local differential 
behavior of the velocity field is then accounted for, at each grid cell, by the symmetric and 
antisymmetric components. The directions of the vorticity, the eigenvectors of the shear 
tensor (which represent directions of collapse or expansion) and halo angular momentum 
have been compared. The main findings are threefold: 

1. Halo spin is strongly correlated with the direction of the local vorticity on scales > 4r v j r . 
The strength of the correlation is mass dependent, with the most massive haloes being 
the most aligned with the vorticity. 

2. The direction of the vorticity tends to be perpendicular to the shear eigenvector asso- 
ciated with the largest eigenvalue (ei). In all web environments except voids (namely 
knots, filaments and sheets) ei corresponds to the axis of fastest collapse. In voids ei 
corresponds to the axis of slowest expansion. 

3. The strength of both alignments depends only weakly on the web classification and 
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the general trends are found throughout the simulated volume. 



A key element of the standard model of cosmology is that structure emerges out of 
rotation-free Gaussian primordial perturbations into a homogenous and isotropic universe. 
It follows that the early dynamical phase of halo grow th is charact erized by being a shear, 
curl-free, (potential) flow. Tidal Torque Theory (TTT, IWhitdll984j ) provides the theoretical 
framework to account for how the angular momentum of a halo's Lagrangian mass grows 
at this early stage. However, as the Lagrangian mass of a halo decouples from the Hubble 
expansion, the flow field around th e halo becomes more non-linear and more vortical (see 
Kitaura et al.ll2012l ). iPorciani et al.l ( 120021 ) showed that at about half the halo's turn around 
time, the growth of the halo's angular momentum no longer follows TTT's predictions. By 
analyzing a high resolution DM-only cosmological simulation we have shown here that at 
the present epoch the spin of simulated halos is strongly aligned with the vorticity of the 
large scale velocity field. This suggests that the acquisition of a halo's spin proceeds in 
two phases, before and after roughly half of the halo's turn around time. In the first phase 
the halo is in its linear and quasi-linear regime and the angular momentum grows out of an 
irrotational shear flow. In the second phase as the halo proceeds towards a virial equilibrium, 
it condenses out of a rotational, vortical flow. At that final stage the spin of a halo reflects 
the vorticity of the flow field of its immediate neighborhood. 

The orientation of galaxies and clusters with respect to the LSS has been shown to 
be an important variable in the process of galaxy formation and has repercussions on for 



example, the correlation function (e.g. see Ivan Daalen et al.ll2012l ). Our findings of a strong 



alignment of halo spin with the vorticity, and the tendency of the vorticity to be orthogonal 
to ei regardless of the web classification, can thus be directly compared with observational 
studies. 



For example, iNavarro et al.l (120041) reported that nearby disc galaxie s (including the 
Milky Way) tend to spin within the super-galactic plane. On greater scales, iBinggelil (119821 ) 
showed that the orientations of clusters of galaxies are anisotropic with respe ct to the cluster 
distrib ution: prolate clusters are roughly parallel. Using the 2MASS survey iLee fc Erdogdu 
( 120071 ) reconstructed the observed shear field finding an alignment between galaxy sp i ns and 
the intermediat e axis. Alignments have been examined in the SDSS by iPaz et al.l (120111 ) 
among others: iTempel et al.l ( 120121 ) found that the spin of spiral galaxies is aligned with 
the filament a xis (i.e. perpendicul ar to ei). The orientation of spin with respect to voids is 
controversial: iTrujillo et al.l ( 120061 ) reported that disc galaxies l ocated in shells s urrou nding 
voids, rotate preferentially in the shell. Thi s was disput e d by ISlosar fc White! (120091 ) who 
found a distribution consistent with random. IVarela et al.l ( 120121 ) find that within specifically 
chosen shells, discs are aligned perpendicular to the void's radial direction. In our web 
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classification scheme voids are bounded by sheets, for which the plane of the sheet is defined 
by ex- Since the vorticity tends to lie within the plane of the sheet , it follows tha t halo spins 
are embedded within sheets - an alignment similar that found by lTrujillo et al.l (120061 ). 



2012 



Much theoretical work (e.g. 


Altav et al. 


2006: 


Patiri et al. 


2006; 


Araeon-Calvo et al. 




2007, 


Brunino et al.ll2007; 


Wang et al. 


2011 


Libeskind et al. 


2012b 


: Codis et al.l2012: 


Trowland et al. 



and references therein) has been dedicated to predicting, understanding and quantify- 
ing these alignments. Although these studies have access to the full six- dimensional phase 
space information, spin alignments with respect to the LSS, but never with the vorticity, 
have until now been investigated. Understanding how the Milky Way halo spins is of direct 



galaxy population ( 


Knebe et al. 


2004; 


Libeskind et al. 


2005. 


2009, 


2011; 


Dcason et al. 


2011 


Pawlowski et al. 


2012). 



The vorticity studied here reflects the behavior of the flow field on scales greater than 
~ 4r vir which for Milky Way type haloes is approximately the Mpc s cale. Just a hand ful of 
dwarf galaxies near the Milky Way have reliable proper motions (see lMetz et al.ll2008l . for a 
description of these): it is thus still too difficult to reliably measure the spin of the Milky 
Way's halo and the vorticity of our cosmic environment. However, future telescopes such as 
LSST and SkyMapper will be able to measure proper motions of dwarfs on these scales, thus 
probing the external flow field. With respect to external galaxies, missions such as GAIA and 
EUCLID may be able to map the large scale orientation of the vorticity in larger volumes. 
This can be done by mapping the orientation of disc galaxies and using it to model the 
orientation of the vortic ity. The well known misa li gnment between s i mulated galactic disks 



and their parent halos ( Bailin fc Steinmetzl 12005 



Bailin et al. 2005; Libeskind et al. 2007 



Bett et al.l 120101 ) complicates the suggested task. Yet, the feasibility of such an attempt 
needs to be tested and gauged against cosmological simulations. 
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Fig. 1. — Left: The probability distribution P(\ cos/i|) as a function of the angle formed, 
| cos/i| between the vorticity and the three principle axes the shear tensor. The distribution 
between oj and ei, e 2 , and e 3 is shown in black, blue and red. The vorticity displays a 
clear perpendicular alignment with ei and parallel alignment with both e 2 and e 3 . The 
average offset between these distributions and a uniform one is ~ 6.5 times the Poissonian 
error of a uniform distribution. Right: The probability distribution between the vorticity 
and the halo spin axis Jhaio- Haloes are binned by mass (10 12 < M/h~ 1 M Q < 10 13 red; 
10 11 < M/h~ l M & < 10 12 blue; and 10 10 < Mjhr x M & < 10 11 black). Poisson error bars for a 
uniform distribution with the same number of haloes in each mass bin are shown along with 
the average offset between the alignment and a random distribution in units of the Poisson 
error (i.e. the statistical significance). None of the alignments are consistent with uniform. 
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Fig. 2. — The orthonormal eigensystem of the shear can be used to define a coordinate 
system. By normalizing the vorticity vectors to unity, their orientation with respect to 
this coordinate system can be shown to be non-uniform. One eighth of the unit sphere 
is shown in an equal area projection here, with contours and colors denoting the (surface) 
density of vorticities. ei points towards the north pole (i.e. cos(#) = 1) while e 2 and 
e 3 are longitudes of and tt/2, respectively. A uniform distribution is represented by a 
monochromatic projection. 
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Fig. 3. — A 250/?r 1 kpc slice through the simulation showing the unit vorticity vector (u) 
(blue arrows) and the principle axis (ei) of the shear tensor. Regions where ei corresponds 
to collapse (expansion) are shown in red (magenta). Since ei is an axis and not a vector, 
arrowheads are not plotted. The density is shown as the grayscale and the white contour 
denotes overdensity with respect to the mean. The length of the ei and u vectors are 
normalized to unity such that a short vector implies perpendicularity with respect to the 
plane. 



